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Abstract 

Spatio-temporal models are widely used for inference in statistics and many applied 
areas. In such contexts interests are often in the fractal nature of the sample surfaces and 
in the rate of change of the spatial surface at a given location in a given direction. In 
this paper we apply the theory of Yaglom (1957) to construct a large class of space-time 
Gaussian models with stationary increments, establish bounds on the prediction errors 
and determine the smoothness properties and fractal properties of this class of Gaussian 
models. Our results can be applied directly to analyze the stationary space-time models 
introduced by Cressie and Huang (1999), Gneiting (2002) and Stein (2005), respectively. 
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1 Introduction 

Spatio-temporal models are widely used for inference in statistics and many applied areas such 
as meteorology, climatology, geophysical science, agricultural sciences, environmental sciences, 
epidemiology, hydrology. Such models presume, on x R, where d is the spatial dimension, a 
collection of random variables X{x,t) at location x and time t. The family {X{x,t) : {x,t) G 
X M} is referred to as a spatio-temporal random field or a space-time model. Many authors 
have constructed various stationary space-time models and the topic has been under rapid 
development in recent years. See, for example, Jones and Zhang (1997), Cressie and Huang 
(1999), de laco, Myers and Posa (2001, 2002, 2003), Gneiting (2002), Gneiting, et al. (2009), 
Kolovos, et al. (2004), Kyriakidis and Journel (1999), Ma (2003a, 2003b, 2004, 2005a, 2005b, 
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2007, 2008), Stein (2005) and their combined references for further information on constructions 
of space-time models and their apphcations. 

There has also been increasing demand for non-stationary space-time models. For example, 
in the analysis of spatio-temporal data of environmental studies, sometimes there is little reason 
to expect stationarity under the spatial covariance structures, and it is more advantageous to 
have space-time models whose variability changes with location and/or time. Henceforth, the 
construction of nonstationary space-time models has become an attractive topic and several 
approaches have been developed recently. These include to deform the coordinates of an 
isotropic and stationary random field to obtain a rich class of nonstationary random fields [see 
Schmidt and O'Hagan (2003), Andcres and Stein (2008)], or to use convolution-based methods 
[cf. Higdon, Swall and Kern (1999), Higdon (2002), Paciorek and Schervish (2006), Calder and 
Cressie (2007)] or spectral methods [Fuentes (2002, 2005)]. 

In this paper, we apply the theory of Yaglom (1957) to construct a class of space-time Gaus- 
sian models with stationary increments and study their statistical and geometric properties. 
The main feature of this class of space-time models is that they are anisotropic in time and 
space, and may have different smoothness and geometric properties along different directions. 
Such properties make them potentially useful as stochastic models in various areas. By apply- 
ing tools from Gaussian random fields, fractal geometry and Fourier analysis, we derive upper 
and lower bounds for the prediction errors, establish criteria for the mean-square and sample 
path differentiability and determine the Hausdorff dimensions of the sample surfaces, all in 
terms of the parameters of the models explicitly. Our main results show that the statistical 
and geometric properties of the Gaussian random fields in this paper are very different from 
those obtained by deformation from any isotropic random field. It is also worth to mention that 
the method in this paper may be applied to analyze more general Gaussian intrinsic random 
functions, convolution-based space-time Gaussian models [Higdon (2002), Calder and Cressie 
(2007)] and the spatial processes in Fuentes (2002, 2005). 

The rest of this paper is organized as follows. In Section 2 we construct a class of space-time 
Gaussian models with stationary increments by applying the theory of Yaglom (1957). Then we 
establish upper and lower bounds for the prediction errors of this class of models in Section 3. 
In Section 4 wc consider smoothness properties of the models and establish explicit criteria for 
the existence of mean square directional derivatives, mean square differentiability and sample 
path continuity of partial derivatives. In Section 5 we look into the fractal properties of these 
models and determine the Hausdorff dimensions of the range, graph and level sets. In Section 
6, we apply the main results to some stationary space-time models, such as those constructed 
by Cressie and Huang (1999), Gneiting (2002) and Stein (2005). Finally, in Section 7, we 
provide proofs of our results. 

We end the Introduction with some notation. Throughout this paper, instead of using 
space-time parameter space M"^ x M, we take the parameter space as or = [0, oo)^. We 
use I • I to denote the Euclidean norm in M . The inner product in is denoted by (•, •). A 
typical parameter, t G M.^ is written as t = (ti, . . . ,tAr). For any s,t £ M.^ such that Sj < tj 
(j = 1, . . . ,N), [s,t] = Ylj=i [sj,tj] is called a closed interval (or a rectangle). 

We will use c, ci , C2 , . . . , to denote unspecified positive and finite constants which may not 
be the same in each occurrence. 
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2 Anisotropic Gaussian models with stationary increments 



We consider a special class of intrinsic random functions; namely, space-time models with 
stationary increments. We will further restrict ourselves to Gaussian random fields for which 
powerful general Gaussian principles can be applied. Many of the results in this paper can 
be extended to non-Gaussian space-time models (such as stable or more general infinitely 
divisible random fields), but their proofs require different methods and go beyond the scope of 
this paper. One can find some information for stable random fields in Xiao (2008). 

Throughout this paper, X = {X[t),t G M^} is a real-valued, centered Gaussian random 
field with -'^(O) = 0. We assume that X has stationary increments and continuous covariance 
function C{s,t) = ¥.[X{s)X(t)]. According to Yaglom (1957), C{s,t) can be represented as 

C{s,t)= [ (e*<^'^> - l)(e-*<*'^> - l)F((iA) + (s,Qt), (2.1) 

where Q is an x non-negative definite matrix and F(dX) is a nonnegative symmetric 
measure on \ {0} satisfying 

/ TTTx\2F(dX)<^. (2.2) 
Jrn 1 + |A|^ 

In analogy to the stationary case, the measure F is called the spectral measure of X. If F is 
absolutely continuous with respect to the Lebesgue measure in R^, its density / will be called 
the spectral density of X. 

It follows from p.ip that X has the following stochastic integral representation: 

X{t)= [ (e'<*'^> -l)H^((iA) + (y,i), (2.3) 

where Xi = X2 means the processes Xi and X2 have the same finite dimensional distributions, 
Y is an A-dimensional Gaussian random vector with mean and covariance matrix Q, W{dX) 
is a centered complex- valued Gaussian random measure which is independent of Y and satisfies 

e(w{A)'W{B)'^ = F{AriB) and W{-A) =W{Aj 

for all Borel sets A,BC M.^ . The spectral measure F is called the control measure of W. 
Since the linear term {Y, t) in (|2.3p will not have any effect on the problems considered in this 
paper, we will from now on assume y = 0. This is equivalent to assuming Q = in ()2.ip . 
Consequently, we have 

v{h) =E{X{t + h) - X{t)f = 2 [ (1 -cos(/i,A))F(dA). (2.4) 

It is important to note that v{h), called variogram in spatial statistics, is a negative definite 
function in the sense of I. J. Schoenberg, which is determined by the spectral measure F. See 
Berg and Forst (1975) for more information on negative definite functions. 

The above shows that various centered Gaussian random fields with stationary increments 
can be constructed by choosing appropriate spectral measures F. For the well known fractional 
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Brownian motion = {B^{t),t € M^} of Hurst index H € (0, 1), its spectral measure has 
a density function 

where c{H, N) > is a normahzing constant such that v{h) = Ih]"^^ . Since v{h) depends on \h\ 
only, is isotropic. Other examples of isotropic Gaussian fields with stationary increments 
can be found in Xiao (2007). We also remark that all centered stationary Gaussian random 
fields can be treated using the above framework. In fact, if y = {Y{t),t € M^} is a centered 
stationary Gaussian random field, it can be represented as Y{t) = f^^ e*<*'^> W{dX). Thus the 
random field X defined by 



X{t) = Y{t)-Y{0)= / (e^<*'^> - 1) W^(dA), VtG 



is Gaussian with stationary increments and X{0) = 0. Note that the spectral measure of X 
in the sense of (12. 4p is the same as the spectral measure [in the ordinary sense] of the stationary 
random field Y. 

In the following, we propose and investigate a class of centered, anisotropic Gaussian ran- 
dom fields with stationary increments, whose spectral measures are absolutely continuous with 
respect to the Lebesgue measure in R^. More precisely, we assume that the spectral measure 
F oi X = {X{t),t G M^} is absolutely continuous with density function /(A) which satisfies 
(j2.2p and the following condition: 

(C) There exist positive constants ci, C2, C3, 7 and • • • , Pn) G (0, 00)^ such that 

^ 1 

7>E^ (2.5) 
j=i J 



and 



< /(A) < 7— xr^^— VAGM^with |A|>C3. (2.6) 



The following proposition shows that (|2.5p ensures / is a legitimate spectral density func- 
tion. 

Proposition 2.1 Assume that /(A) is a non-negative measurable function defined on M^. // 

|A|2/(A)dA < 00 

|A|<1 

and \2. 6]) holds, then /(A) is a legitimate spectral density if and only if the parameters 7 and 
Pj for j = !,■■■ ,N satisfy (Ej 



Some remarks about condition (j2.6p are in the following. 
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Remark 2.2 



There is an important connection between the random field models that satisfy Condition 
(C) and those considered in Xiao (2009). For j = 1, - ■ ■ ,N, let 



and let Q = J2f=i TT' Then (j2.6p can be rewritten as 



^Wtq ^ /(^) ^ ^^N ,?,H,^2+Q - V A G with |A| > C3, (2. 



where the positive and finite constants C4 and C5 depend on N, ci,C2, /3j and 7 only. 
To verify this claim, we will make use of the following elementary fact: For any positive 
numbers and q, there exist positive and finite constants C4 and C5 such that 

N a N / N 



for all non-negative numbers ai, . . . , oat. Note that 



■ 2 W 



and i(7 - X^^Ii ;J-)(2 + Q)=1- We see that ([IJ]) and are equivalent. 

In turns out that the expression (j2.8p is essential in this paper and will be used frequently. 
For simplicity of notation, from now on we take C3 = 1. 

It is also possible to consider Gaussian random fields with stationary increments whose 
spectral measures are not absolutely continuous. Some examples of such covariance 
space-time models can be found in Cressie and Huang (1999), Gneiting (2002), Ma 
(2003a, 2003b). Since the mathematical tools for studying such random fields are quite 
different, we will deal with them systematically in a subsequent paper. 

Nonstationary Gaussian random fields can be constructed through deformation of an 
isotropic Gaussian random field. Refer to Anderes and Stein (2008) for more details. One 
of the advantages of deformation is to closely connect a nonstationary and/or anisotropic 
random field to a stationary and isotropic one for which the existing statistical techniques 
are available. However, there is also a disadvantage [from the point of view of flexibility] 
associated with deformation. Let X{t) = Z{g^^{t)), where {Z{t),t G M^} is an isotropic 
Gaussian model and 5 is a smooth bijection of R^. Since the function g is bi-Lipschitz 
on compact intervals, the fractal dimensional properties of X are the same as those of 
Z. Hence deformation of isotropic Gaussian models will not generate anisotropic random 
fields with rich geometric structures as shown by the models introduced in this paper. 
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3 Prediction error of anisotropic Gaussian models 



Suppose we observe an anisotropic Gaussian random field X on at t^, . . . and wish to 
predict X{u), for u € M^. Tlien the inference about X{u) will be based upon the conditional 
distribution of X{u) given the observed values of X{t^), . . . , Xif^). Refer to Stein (1999, 
Section 1.2) for the closed form of this conditional distribution. A statistical analysis typically 
aims at the optimal linear predictor of this unobserved X{u), known as simple kriging. The 
simple kriging predictor of X{u) is 

X*{u) = c{uf-E-^Z, (3.1) 

where Z = (X(ti), . . . , X(r))^, c{uf = Cov{X{u),Z} and E = Cov(Z,Z^). The form (l3J]) 
minimizes the mean square prediction error, which then is given as Var(X(n))-c(n)^S-^c(n). 
Since X is Gaussian, the simple kriging is the conditional expectation of X(u) given Z, and 
the mean square prediction error is the conditional variance of X(u) given Z. 

The main result of this section is Theorem 13. II below, which gives lower and upper bounds 
for the mean square prediction error for Gaussian random fields with stationary increments 
which satisfy Condition (C). It shows that, similar to stationary Gaussian field models [cf. Stein 
(1999)], the prediction error of the models in this paper only depends on the high frequency 
behavior of the spectral density of X. 

Theorem 3.1 Let X = {X{t),t € M^} be a centered Gaussian random field valued in M with 
stationary increments and spectral density /(A) satisfying ^2. 6\) . Then there exist constants 
cq > and cj > 0, such that for all integers n > 1 and all u,t^, - ■ ■ , E M^, 

N N 
C6 min < Var(X(tx)lX(ti),... ,X(t")) <C7 min - tj\) , (3-2) 

where Hj is given in \2. 7| ), t'^ = and aj : — is defined by 

r^^^ if 0< Hj < 1, 

C7,(r) = { r^\logr\ if Hj = 1, (3.3) 



„2 



if H, > 1. 



If Hj < 1, for j = 1, • • • ,N, then the two bounds in ()3.2p match. When there is some 
Hj > 1, that means, the random field X{t) is smoother in the j-th direction [see Gorollaries 
14.31 14.71 and Theorem 14.81 below], then the upper and lower bounds are not the same any more. 
This suggests that the prediction error may become bigger as X{t) gets smoother in some 
directions. 

The proof of Theorem 13. H as well as those of Theorems 14.91 15.11 and 15.21 reply partially on 
the following lemma, which provides upper and lower bounds for the variogram of the model. 

Lemma 3.2 Let X = {X{t),t G M^} be a centered Gaussian random field valued in M. with 
stationary increments and spectral density /(A) satisfying h2. Then there exist constants 
cg > and cg > such that for s,t ^ R^, 

N N 

C8^f^j(|s, -til) <E{X{s) - X{t))^ <cgY,aj{\sj -tj\), (3.4) 
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where the function cTj is defined in \3. 3)) . 



The upper bound in (j3.4p implies that X has a version whose sample functions are almost 
surely continuous. Throughout this paper, without loss of generality, we will assume that the 
sample function 1 1— > X{t) is almost surely continuous. 

4 Smoothness properties of anisotropic Gaussian models 

Regularity properties of sample path of random fields are of fundamental importance in proba- 
bility and statistics. Many authors have studied mean square and sample path continuity and 
differentiability of Gaussian processes and fields. See Cramer and Leadbetter (1967), Alder 
(1981), Stein (1999), Banerjee and Gelfand (2003), Adler and Taylor (2007). In this section we 
provide explicit criteria for the mean square and sample path differentiability, for the models 
introduced in Section 2. 

4.1 Distributional properties of mean square partial derivatives 

Banerjee and Gelfand (2003) studied the smoothness properties of stationary random fields and 
some non-stationary relatives through directional derivative processes and their distributional 
properties. To apply their method to random fields with stationary increments, let us first 
recall the definition of mean square directional derivatives. 

Definition 4.1 Let u € be a unit vector. A second order random field {X{t),t G M^} has 
mean square directional derivative X'^{t) at t ^ in the direction u if, as h ^ 0, 

Xit + hu)-X{t) 

-^u,h{'^) = 

converges to X'^{t) in the L2 sense. In this case, we write X'^{t) = l.i.m.ft^o^u,h(i)- 

Let 61,62,- •• ,6Ar be an orthonormal basis for M^. If ti = ej, then X'^.{t) is the mean 
square partial derivative in the j-th. direction defined in Adler (1981), which will simply be 
written as Xj{t). We will also write X^.^^^t) as Xjh{t). 

For any second order, centered random field {X{t),t G M^}, similar to Theorem 2.2.2 in 
Adler (1981), one can easily establish a criterion in terms of the covariance function C{s,t) = 
E[X(t)X(s)] for the existence of mean square directional derivative X'^{t). Banerjee and 
Gelfand (2003) further showed that the covariance function of X'^{t) is given by 

Ku{s,t) = lim limE[Xu,h{t)Xu,kis)] 

C(t + hu,s + ku)-C(t + hu,s)-C(t,s + ku) + C(t,s) 
= lim lim — . 

h^O k^O hk 

Extending their argument, one obtains the following theorem for Gaussian random fields with 
stationary increments. 
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Theorem 4.2 Let X = {X{t),t € M^} he a centered Gaussian random field valued in M with 
stationary increments, then the mean square partial derivative Xj{t) exists for all t G if 
and only if the limit 

j.^ vjhej) + vjkej) -v{{h- k)ej) 
h,k^o hk 

exists. Moreover, this later condition is equivalent to v{t) has second-order partial derivatives 
at in the j-th direction. 

As a consequence, we obtain an explicit criterion for the existence of mean square partial 
derivatives of Gaussian random fields in Section 2. 

Corollary 4.3 Let X = {X{t),t G M^} be a centered Gaussian random field valued in M 
with stationary increments and spectral density /(A) satisfying Condition (C). Then for every 
j = 1, - ■ ■ ,N, the mean square partial derivative ^jit) exists if and only if 

' ' > 2, (4.2) 



or equivalently Hj > 1 [cf. ^2. 7p ] . 

Assume condition (j4.ip of Theorem 14.21 holds so that the mean square partial derivative 



exists for all t G M . We now consider the distributional properties of the random field 
{Xj(t),tGM^}. 

Since = for ah t G M^, we have E(Xj- = and = 0. Let 

Cj^\s,t) and Cj{s,t) denote the covariance functions of the random fields {Xj^h{t),t G M^} 
and {Xj{t),t G M^}, respectively. Let A = s — t, we immediately have 

•' 2h^ 
and Yar{Xj^h{t)) = v{hej) /h^ , which only depends on the scalar h. 

Theorem 4.4 Let X = {X{t),t G M^} he a centered Gaussian random field valued in M with 
stationary increments. Suppose that all second-order partial derivatives of the variogram v{t) 
exist. Then the covariance function of X'- (t) is given by 

C,{s,t) = y;{s-t), (4.4) 

where Vj (t) is the second-order partial derivative of v at t in the j-th direction. In particular, 
{Xj{t),t G M^} is a stationary Gaussian random field. 

Proof The desired result follows from (14.31). □ 
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It is also useful to determine the covariance of X{s) and X'j{t) for all s,t G M^. Since 

Coy{Xis),Xj,h{t)) = ^ {vit + hej) - vit) + v{A) - viA - hej)}, 
where A = s — t, we obtain 

Cov{Xit),Xj,h{t)) = ^ {vit + hej) - vit) - vihej)] 



and 



Cov(X(s),Xj(i)) = lim ^ {vit + hej) - vit) + viA) - viA - hej)} 



1 



v'^it) + v'^iA)), 



(4.5) 



where v'^it) is the partial derivative of v at t in the j-th direction. 

In particular, Cov(X(t), Xj- (t)) = v'jit)/2, which is different from the stationary case. 
Recall that if Yit) is a stationary Gaussian field with mean square partial derivative Yjit), 
then Yit) and Y^it) are uncorrected. That means, the level of the stationary random field at a 
particular location is uncorrclatcd with the partial derivative in any direction at that location. 
However, this is not always true for nonstationary random fields. 

Next we consider the bivariate process 

It can be verified that this process has mean and cross-covariance matrix 

yj,hi^,t) 

/ vis) + vit) - viA) vit + hej) - vit) + viA) - viA - hej) \ 

^ vis + hej) -vis) + viA) - viA + hej) viA + hej) + v^^ - hcj) - 2viA) 
\ 2h 2^2 / 

Because Yj^\t) is obtained by linear transformation of the above is a valid cross- 

covariance matrix in R^. Since this is true for every h, letting h ^ we see that 



Vjis,t) 



^{vis)+vit)-viA)} ^{v'^it)+v'^iA)} ^ 
l{v'^is)-v'jiA)} l.;(A) ^ 



is a valid cross-covariance matrix in R^. In fact, Vj is the cross-covariance matrix for the 
bivariate process 

Xit) 



Yjit) 



x'Jt) 
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4.2 Criterion for mean square differentiability 

Benerjee and Gelfand (2003) pointed out that the existence of all mean square directional 
derivatives of a random field X does not even guarantee mean square continuity of X, and 
they introduced a notion of mean square differentiability which has analogous properties of 
total differentiability of a function in in the non-stochastic setting. We first recall their 
definition. 

Definition 4.5 A random field {X{t),t G M^} is mean square differentiable att ^ if there 
exists a (random) vector Vx{t) € such that for all scalar h > 0, all vectors u S Sn = {t ^ 
: |t| = 1} 

X{t + hu) = X{t) + hu^Vx{t) + r{t, hu), (4.6) 
where r(t, hu)/h in the L2 sense as h ^ 0. 

In other words, for all vectors u E Sn, it is required that 

h^o \ h J 

It can be seen that if X is mean square differentiable at t, then for all unit vectors u € Sn 

l.i.m. X{t + hu) - X{t) 
^uK^) = ^-"^ 

l.i.m. hu^Vxit) + r{t,hu) 

= h^O ; 

h 

= u^Vx{t). 

Hence it is necessary that Vx(i) = {X[{t), . . . ,X'^{t)). 

The next theorem provides a sufficient condition for a Gaussian random field with stationary 
increments to be mean square differentiable. 

Theorem 4.6 Let X = {X{t),t € M^} he a centered Gaussian random field valued in M with 
stationary increments. If all the second-order partial and mixed derivatives of the variogram 
v{t) exist and are continuous, then X is mean square differentiable at every t G M^. 

As a consequence of Theorem 14.61 we obtain 

Corollary 4.7 Let X = {X{t),t G M^} be a centered Gaussian random field valued in M with 
stationary increments and spectral density /(A) satisfying Condition (C). If 

N 

I \ 

> 2 for every j = 1, . . . ,N, (4.8) 



then X is mean square differentiable at every t G M.^ . 
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4.3 Criterion for sample path differentiability 



For many theoretical and applied purposes, one often needs to work with random fields with 
smooth sample paths. Refer to Adler (1981), Adler and Taylor (2007) and the reference therein 
for more information. Since in general mean square differentiability does not imply almost sure 
sample path differentiability, it is of interest to provide convenient criteria for the latter. 

For Gaussian random fields considered in this paper, it turns out that under the same 
condition as Corollarv 14.31 the partial derivatives of X are almost surely continuous. 

Theorem 4.8 Let X = {X{t),t G M^} be a separable and centered Gaussian random field 
with values in M. We assume that X has stationary increments and satisfies Condition (C). 
(i). If 

' ' > 2 {i.e.,Hj > 1), (4.9) 



for some j G {1, • • • , A^}, then X has a version X with continuous sample functions such that 
its jth partial derivative Xj{t) is continuous almost surely. 

(ii). // i4-^ holds for all j G {1, • • • ,N}, then X has a version X which is continuously 
differentiable in the following sense: with probability 1, 

Xit + hu) - X(t) -hu^V ^(t) J.J 

lim — ^ ^ = for all u (£ Sn andteR^. (4.10) 

h^O h 

If condition (j4.9p does not hold for some j G {1, • • • , N}, then X{t) does not have mean 
square partial derivatives along those directions and X{t) is usually a random fractal. In this 
case, it is of interest to characterize the asymptotic behavior of X{t) by its local and uniform 
moduli of continuity. 

These problems for anisotropic Gaussian random fields have been considered in Xiao (2009) 
and the methods there are applicable to X with little modification. For completeness, we state 
the following result which can be proved by using Lemma 13.21 and general Gaussian methods. 
We omit its proof. 

Theorem 4.9 Let X = {X{t),t G M^} be as in Theorem \4-8\ Then for every compact 
interval I C M^, there exists a positive and finite constant cio, depending only on I and 
Hj, (j = 1, . . . , N) such that 

supie^ g[o,e] |X(t + s) -X(t)| 
limsup = < cio, (4.11) 

where ip{e) = Ylf=i ^jil^jl) f^'^ ^ = (^i; ■ ■ ■ i^n) G I^^- 



5 Fractal properties of anisotropic Gaussian models 

The variations of soil, landform and geology are usually highly non-regular in form and can be 
better approximated by a stochastic fractal. Hausdorff dimension has been extensively used 
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in describing fractals. We refer to Kahane (1985) or Falconer (1990) for their definitions and 
properties. 

Let X = {X(t),t € M^} be a real- valued, centered Gaussian random field. For any integer 
p>l,we define an (A^,p)-Gaussian random field X = {X(t),t G R^} by 

X(t) = (Xi(t),...,Xp(t)), ^GM^ (5.1) 

where Xi, . . . , Xp are independent copies of X. 

In this section, under more general conditions on X than those in Sections 2-4, we study the 
Hausdorff dimensions of the range X([0, 1]^) = {X(t) : t e [0, 1]^}, the graph GrX([0, 1]^) = 
{{t,X{t)) : t G [0, 1]^} and the level set X-i(x) = {t G : X{t) = x} {x e W). The results 
in this section can be applied to wide classes of Gaussian spatial or space-time models (with 
or without stationary increments). 

First, let's consider fractional Brownian motion = {B^{t),t G M^} valued in MP with 
Hurst index H G (0, 1). B^{t) is a special example of our model which, however, has isotropic 
spectral density. It is known [cf. Kahane (1985)] that 



a.s. 



dimGr5^([0,l]^) = min|iV+ (1 -i/)p, ^| 
Especially, when p = 1, 

dimGrB^([0,l]^) = A^ + l-iJ a.s. 

Also, 

dim {B")-^{x) =N-H, a.s. 

The fractal properties of fractional Brownian motion have been applied by many statisticians 
to estimate the Hurst index H and it is sufficient to choose p = 1. Refer to Hall and Wood 
(1993), Constantine and Hall (1994), Kent and Wood (1997), Davis and Hall (1999), Chan and 
Wood (2000, 2004), Zhu and Stein (2002). 

Let {Hi, . . . ,Hn) € (0, 1]^ be a constant vector. Without loss of generality, we assume 
that they are ordered as 

< < :ff2 < • • • < ^jv < 1. (5.2) 

We assume the following conditions. 

(Dl). There exist positive constants Sq, cn > 1 such that for all s,t e [0, 1]-^ with \s—t\ < So 

TV _ ^ _ 

i=i i=i 

(D2). For any constant e G (0, 1), there exists a positive constant ci2 such that for all u,t 
G [e, 1]^, we have 

Var {X{u) I X{t)) > cu E - tjf^'- (5-4) 

The following theorems determine the Hausdorff dimensions of range, graph and level sets 
of X. Because of anisotropy, these results are significantly different from the aforementioned 
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results for fractional Brownian motion or other isotropic random fields [cf. Xiao (2007)]. Even 
though Theorems 15.11 and 15.21 below are similar to Theorems 6.1 and 7.1 in Xiao (2009), they 
have wider applicability. In particular, they can be applied to a random field X which may be 
smooth in certain (or all) directions. 

Theorem 5.1 Let X = {X(t),t € M^} be an {N,p)- Gaussian random field defined by i5. 
If the coordinate process X satisfies Condition (Dl), then, with probability 1, 



r ^11 

dimX([0,l]^) =min p; 



and 

dimGrX([0, 1]^) = ^min^ {^^+^"^ + (1- ^''^^ E ^ }' 
where Y.]=i ^ ■= 0- 



(5.5) 



(5.6) 



Proof The right inequality in (15. 3p and Theorem 4.9 show that X(i) satisfies a uniform 
Holder condition on [0, 1]^ which, in turn, imply the desired upper bounds in ()5.5p and (j5.6p . 

The lower bounds for dimX([0, 1]^) and dimGrX([0, 1]^) can be derived from the left 
inequality in (jS.Sp and a capacity argument. See the proof of Theorem 6.1 in Xiao (2009) for 
details. □ 

For the level sets of X, we have 

Theorem 5.2 Let X = {X(t),t € M^} be an {N,p)- Gaussian random field defined by i5.1\) . 
If the coordinate process X satisfies Conditions (Dl) and (D2), then the following statements 
hold: 

(i) // Y.f=i < ^/^en for every x e W\{0}, X-i(x) =0 a.s. 

(ii) // X^^i =- > p, then for any x € W, with positive probability 

dimX~^(x) = min I =^ + N — k — Hh 




Proof The results (i) and (ii) follow from the proof of Theorem 7.1 in Xiao (2009). □ 

Let X = {X{t),t € R^} be a centered Gaussian random field valued in M with stationary 
increments and spectral density /(A) satisfying ()2.6I1 . Let Hi, . . . , Hjy be defined by ()2.7p . 
Then, by Theorem 13.11 and Lemma 13.21 we see that X satisfies (Dl) for all Hj G (0, 1 A Hj) 
(1 < j < N). It also satisfies Condition (D2) with Hj = Hj provided Hj < 1 for all j = 
1, • • • ,N. Hence one can apply Theorems 15.11 and 15.21 to derive the following result. 

Corollary 5.3 Let X = {X(t),t S M^} be a centered Gaussian random field valued in MP 
defined by i5. We assume that its coordinate process X has stationary increments with 
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spectral density /(A) satisfying 112. 6\) and Hj (j = 1, • • • , N) defined by \2. 7| j are ordered as 
Hi<H2<---< Hn. We have 

(i) . With probability 1, (ESjj and l^j hold with Hj = IA Hj (I <j < N). 

(ii) . //, in addition, Hj < 1 for all j = 1, • • • and J2jLiW' ^ P' ^^^^ 7| j /loZds wi^/i 
positive probability. 

We believe that the above fractal properties can also be useful for estimating the parameters 
Hi, . . . , Hn of our model. However this will be more subtle than the isotropic case, where only 
the single parameter is involved, for the following two reasons. First, if a parameter Hj > 1, 
then the sample function X{t) is smooth in the jth direction and the Hausdorff dimensions of 
X has nothing to do with Hj. In other words, based on fractal dimensions, a parameter Hj 
can be explicitly estimated only when Hj < 1. 

Secondly, if we let p = 1, then (15. 6p gives dimGrX([0, 1]^) = N + 1 - Hi, which does 
not give any information about the other parameters H2, . . . , H^. This suggests that, in order 
to estimate all the parameters of an anisotropic random field model, one has to work with a 
multivariate random field X as defined by (|5.1|) . 



6 Applications to some stationary space-time models 

The above results can be applied to the stationary space-time Gaussian fields constructed by 
Cressie and Huang (1999), Gneiting (2002), de laco, Myers, and Posa (2002), Ma (2003a, 
2003b) and Stein (2005). 



6.1 Stationary covariance models 

Extending the results of Cressie and Huang (1999), Gneiting (2002) showed that, for {x,t) G 

■md 



X 



Cj2 / C|X|27 



is a stationary space-time covariance function, where o" > 0, a > 0, c > 0, a G (0, 1], /3 G (0, 1] 
and 7 G (0, 1] are constants. It can be verified that the corresponding spectral measure is 
continuous in space x and discrete in time t. See Ma (2003a, 2003b) for more examples of 
stationary covariance models. 

In the following, we verify that the sample functions of these space-time models are fractals. 
We will check Conditions (Dl) and (D2) first, and then obtain the corresponding Hausdorff 
dimension results from Theorems 15.11 and 15.21 



Proposition 6.1 Let X = {X{x, t), {x, t) G R'^xM} be a centered stationary Gaussian random 
field in M with covariance function as 116. Then for any M > 0, there exist constants C13 > 
and ci4 > such that 

ci3{\x - + \t- < E{X{x,t) - X{y, s)f < cm{\x - y\^^ + \t - sp") (6.2) 

and 

Yar{X{x,t)\X{y,s)) > C13 (jx - + |t - s^") (6.3) 
for all {x,t) and {y,s) G [-M,M]'^+^. 
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Proposition 6.2 Let X = {X{x, t), [x, t) G M'^ xM} be a centered stationary Gaussian random 
field in M with covariance function as \6. and let X he its associated {N,p)-random field 
defined by Ii5.1\) . Then, with probability 1, 



dimX([0, 1]''+^) = mill <{ p; - + 



d 



7 a 



(6.4) 



And i/ < a < 7 < 1, then 

dimGrX([0,l]'^+i) 
//O < 7 < a < 1, then 

dimGrX([0, 1]"'+^; 



d+l + (l-a)p 
d+l + {l-l)p 



I , d 
a 7 



d+l + {l--f)p 
^ + l + {l-a)p 

a 7 



if i < p < i + ^, 

a — ^ 07' 
if p > 1 + rf. 



if f <P<^ + f, 



(6.5) 



(6.6) 



Proposition 6.3 Let X = {X(x, t), (x, t) € M^xM} be a centered stationary Gaussian random 
field in M with covariance function as ^6.1\) . and let X he its associated {N ^p) -random field. 
(i) When ^ + ^ < then for every x £ M^, X~^(x) = a.s. 



d 

a ' 7 

(ii) When ^ + ^ > P; ^/0<a<7< 1, t/ien for any x G M*', with positive probability 



dimX ^(x) = 
and i/ < 7 < a < 1, then for any x G 

dimX^^(a;) = 



d + 1 — ap 



with positive probability 



d + 1 — ^p 
+ 1 — ap 



da 



(6.7) 



(6.8) 



6.2 Stationary spectral density models 

In Section 6.1, the stationary space-time models are constructed directly by covariance func- 
tions, which are isotropic in the space variable. Stein (2005) showed that stationary covariance 
functions which are anisotropic in space can be constructed by choosing spectral densities of 
the form 



fW={Y.^c,{a, + \X,\r'^ 



where > 0, Cj > 0, aj > and aj G N for j = 1, • • • 

d+l 

y-<2u. 



VA G M'* X M, (6.9) 
d+l are constants such that 



This last condition guarantees / G L^(M'^+^). Clearly /(A) in (j6.9p satisfies (j2.6p with f]j = aj 
and 7 = 2i/. Hence we may apply our results to analyze this class of models, through the 
smoothness properties and the fractal properties. 
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Proposition 6.4 Let X = {X{x, t), {x, t) G M'^ xM} be a centered stationary Gaussian random 
field in M with spectral density as lid. 9\) . 

(i) If 

d+l ^ 

> y — + 

i=i 

then X{x,t) is mean square differentiahle and has a version X{x,t) which is sample path 
differentiable almost surely. 

(ii) X is a fractal [i.e. the sample path of X may have fractional Hausdorff dimension] if and 
only if 

d+i ^ d+i ^ 2 

V— <2z^<y — + — . 

~[ «i ~^ «i mini<£<d+i ae 

The Hausdorff dimensions of various fractals generated by this kind of models can also be 
computed using Corollary 15.31 with Hj = aj{v — X]^=i 207) ' -^i = 1 A //j for j = 1, • • • , d + 1. 
We leave the details to an interested reader. 



7 Proofs 

Proof of Proposition 12.11 

Since (|2.2p is equivalent to J^n (l A |A|^)/(A)dA < 00, and Jj;^|<;^ \X\^ f {\)d\ < 00 is given, 
it is enough for us to show 

d\ 

is equivalent to (|2.5p . 

For this purpose, we appeal to the following fact: Given positive constants (3 and 7, there 
exists a finite constant C15 such that for all a > 0, 

/ ci5 a-(^-^^ if /37 > 1, ^^ -^^ 



lo {a + xf^p I +00 if /37 < 1. 

To verify this, we make a change of variable x = a^y to obtain 

Jo {a + x(^p Jo {l + yl^p- 

Thus (HH) follows. 

First we assume ()2.5p holds. Since | A| > 1 implies that | Aj^ | > for some jo G {1, ... , N}. 
Without loss of generality we assume jo = 1. Then by using (17.1[) (A^ — 1) times we obtain 

dX AT /"°° d\2 • • • dX]\f—i 



dXi 

Vw (lAi|/3i)^~^^=^^ 
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because /3i(7 — Yl^=2 > ^- This proves the sufficiency of (j2.5p . 



To prove the converse, we assume (|2.5|) does not hold. Then there is a unique integer 

AT 



r G {1, . . . , iV} such that Ya=i < ^ ^ ELi 3-- Note that 



By using (j7.ip and integrating dAi-'-dA,-, we see that the last integral is divergent. This 
finishes the proof. □ 

Proof of Lemma 13.21 

For any s,t e M^, denote sq = t, si = {si,t2,--- ,tN), h = {si, S2,t3, ■ ■ ■ ,1^), 
sn-1 = (si, • ■ ■ ; sj\[-i,tj\[) and sn = s. Let h = s — t = {hi, • • • , hj\[). By Jensen's inequality, 
([2^ and ([Zl]) we can write 

N 



E{X{s) - X{t)y < N ^E(X(sfe) - X{sk-i)y 

k=l 

N 

= 2ivV/ (l-cos(/ifcAfc))/(A)dA 



<2ivV/ (l-cos(/ifcAfc))/(A)dA 

fe=l^|A|<l 

^ /■ dA 



(7.2) 



= /l+/2. 



1 -r -'2- 

2 



By using the inequality 1 — cos x < x we have 

N 

|A[2/(A)dA < ci6|s-tp (7.3) 



\tr[ J J\x\<i 



for some positive and finite constant cig. 

To bound the fcth integral in I2, we note that, when |A| > 1, either |Afc| > or there is 

Jo 7^ k such that \Xjg\ > We break the integral according to these two possibilities. 
[1 - cos{hk\k))- ^"^ 



^ o /"°° /'1 \ \\^\ f dXi ■ ■ ■ dXk-idXk+1 ■ ■ ■ dXN 



(7.4) 







+ 4/ {l-cos{hkXk))dXk I ^ dX 
h + h, 
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where dX^j^ denotes integration in Aj {i ^ k^jo). 

By using (jT.ip repeatedly [A^ — 1 times], we obtain 



oo 



1 - cos{hkXk) 



/liAl , /■°° 1 ^^ \ (7-5) 



< cak{\hk\), 



where ak is defined as in p.3[) . 

Similarly, we use (j7.1|) N — 2 times to get 



l4<c J (l - cos{hkXk))dXk j 



dXjo 



< c 



/' (1 - cos(/i,A,))dA, (7-6) 



^0 



<c|/lfc|2. 



Combining (j7.2|) - (j7.6|) yields the upper bound in (j3.4|) . 

Next we prove the lower bound in (j3.4p . By (|2.4p and (j2.8p we have 

E{Xis) - Xit))' > C4 ^^^^^ (1 - cos - t, A)) ^p^, (7.7) 

where p{X) = J2f=i , ^ e 1^^- So, for the lower bound of E{X{s) - it is enough 

to show that for every j = 1, ■ ■ ■ ,N and all h € M^, we have 

/ (1 - cos (/i. A)) >ccT,(|/ij|), (7.8) 

J\\\>1 p[A)'^^ 

where c is a positive constant. 

We only prove (j7.8p for j = 1, and the other cases are similar. Fix h S with \hi\ > 
[otherwise there is nothing to prove] and we make a change of variables 

We consider a subset of the integration region defined by 

D{h) = |y G : \yi\ G [p(/i)^i~\ 1], < 1 and yihi > for 1 < £ < ivj. 
Since p{X) = p{y)/p{h), we have 

dX 



^ — cos ( y^f. hep(h) I yi , 
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By using the inequality 1 — cos x > cx^ for all |x| < A^, where c > is a constant, and the fact 
that hiyi > for all 1 < ^ < A^, we derive that the last integral is at least [up to a constant] 



>p{hf / .hlp{h) ^lyjdy 



> cp{h) «i hi 



^ I- 2 , /"^ r dy2---dyN 

^ yfdyi 



^ n^pyu) -^yiayi ■■■ |f7^Q+2 

1 ioJ^{Ee=i\yi\^') (7.10) 



N-1 



= cai{\hi\). 

This proves (j7.8p and hence Lemma |3.2[ □ 

In order to prove Theorem 13 -H we will make use of the following lemma which implies that 
the prediction error of X is determined by the behavior of the spectral density /(A) at infinity. 



Lemma 7.1 Assume ^2. 0i) is satisfied, then for any fixed constant T > 0, there exists a positive 
and finite constant cn such that for all functions g of the form 



ff(A)=j;a,(e^<*'='^)-l), 

k=l 

where € M and t^ G [—T,T]^, we have 

|<7(A)| <ci7lA|( / \gm^fmdi ' (7.11) 



1/2 



for all A G that satisfy |A| < 1. 

Proof By ()2.6p . we can find positive constants C and rj, such that 

C 

fW>TrT-, VA G with |A| large enough. 

\ \ 

Then the desired result follows from the proof of Lemma 2.2 in Xiao (2007). □ 

Proof of Theorem 13.11 

First, let's prove the upper bound in (j3.2p . By Lemma 13.21 we have 

Yai(x{u)\X{t^),--- ,X(r)) < min E(X(n) - X(t'=))^ 

V / 0<k<n 



N (7.12) 



0<fc<n- 
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In order to prove the lower bound for the conditional variance in ()3.2p . we denote r = 

niin y^^i \uj — t^\^^ . Working in the Hilbert space setting, the conditional variance is just 

0<fc<n ■' 

the square of L^(P)-distance of X{u) from the subspace generated by {X(t^), ■ • • ,X{t'^)}, so 
it is sufficient to prove that for all € M, 1 < A: < n, 



^ k=i ' 



(7.13) 



where ce is a positive constant which may only depend on . . . , Hn and N. 

By using the stochastic integral representation (j2.3p of X, the left hand side of (j7.13p can 
be written as 



E(x{u)-y2akX{t^)) = [ 



J{u,X) 



k=l 



fWdX. 



(7.14) 



Hence, we only need to show 



fc=0 



/(A) d\>c^r\ 



(7.15) 



where = and oq = 1 — X]fc=i '^k- 

We choose a function 6{-) : [0,1] in C°°(]R^) [the space of all infinitely differ- 

entiable functions defined on M^] such that 5(0) = 1 and it vanishes outside the open set 
|t £ : X^jLi < l}- Denote by 6 the Fourier transform of 6. Then one can verify 

that j(-) G C°°(M^) as well and (5(A) decays rapidly as |A| — > cxd. 

Let E be the N x N diagonal matrix with H^^ , ■ ■ ■ ,H^^ on its diagonal and let Sr{t) = 
r-Q 8{r-^t) for ah t G M^. Then the inverse Fourier transformation and a change of variables 
yield 



5^{t) = (27r) 



~N 



-Kt^)8{r'^\) dX. 



(7.16) 



Since min { ^f=i Wj — tj\^^ : < A; < n| > r, we have 6r{u — t'^) = for k = 0,1, ■ ■ ■ ,n. This 
and (j7.16p together imply that 



J 



,A> ^-i(u,X) 



6{r^X) dX 



k=0 

n 



(7.17) 



= (2^)^U,(0)-^afc5,(7x-t'=) 

V fc=o / 

= (2^)^r-«. 

Now we split the integral in (j7.17p over {A : |A| < 1} and {A : |A| > 1} and denote the two 
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integrals by Ii and I2, respectively. It follows from Lemma |7. II that 



h < 



'|A|<1 

< Ci7 

< C18 



J{u,X) 



E 



k=0 



\6{r^X)\dX 



J{u,X) 



E 

fc=0 



A) 



E[x(u)-f;afcX(t'=)] 
. ^ fc=i ^ 





1/2 


'/(A)<iA 


/ |A||5(r^A)|dA 




^|A|<1 


1/2 





(7.18) 



where the last inequality follows from ()7.14p and the boundedness of 6. 

On the other hand, by the Cauchy-Schwarz inequality and (|7.14p . we have 



-'2 ^ 



|A|>1 



J(u,X) 



fc=o 



i(t'=,A> 



fWdX 



1 



|(5(r-^A)pdA 



>i /(r-^A) 



/ " \ 2 

<E(x{u)-Y,akX{t'')] r-' 

^ fc=i ^ 



|A|>1 /(A) 

1 



|5(A)rdA 



(7.19) 



\SiX)\'^dX. 



The last integral is convergent thanks to the fast decay of 6{X). Finally, combining ()7.17p . 
(msll and (I719I) . we get 



(27r)^r-^ < C19 



k=l 



1/2 



Henceforth (17.130 follows, and the theorem was proved because of (17.120 and ()7.13p . 



□ 



Proof of Theorem 14.21 

For t € M^, it is known that Xj^h = '^^^^ converges in L^, as /i ^ 0, if and only if 

Dh,k = {Xit + hej) - X(t)) {X{t + kej) - X{t)) } 

converges to a constant as /i, /c — > 0. However, 

Dh,k = -^\c{t + hej,t + kej) -C{t,t + kej) -C{t + hej,t) + C{t,t)\ 



(7.20) 



So the first part of the theorem is proved. For the second part, it is clear that if v{t) has 
second-order partial derivatives at in the j-th direction then (j4.ip holds [thanks to Taylor's 
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theorem]. On the other hand, if (j4.ip holds, then by taking /i = — s- in ()7.20p we see that 
dv/dtj{0) = 0. This fact, together with ()4.ip . imphes that 

d'^v , 1, v{{k + h)en) - v{ken) 
—^ 0) = hm - hm '-^ 

dtj^ fc^o k h-^Q h 



hm hm 



v{{k + h)ej) — v{kej) + v{hej) 



k^o h^o hk 
exists. This finishes the proof of Theorem 14.21 □ 

Proof of Corollary 14.31 

By Theorem 14.21 it amounts to show that hm Dh^k exists if and only if (14. 2p holds, [i.e., 

h,k—>0 ' 

Pjh - Eili j-) > 2-] It follows from ([7:20]) and that 

0^1^=1 1 - cos {hej, A) - cos {kej,X) + cos {{h - k)ej, A) ^^^.^ 
Jm hk 

To prove the sufficiency of ()4.2p . we note that for each fixed A € M^, 

1 — cos(/iAj) — cos(A;Aj) + cos((/i — A;)Aj) 2 

hm — — A,- 

h,k^O hk ■' 

and by the mean value theorem. 



(7.22) 



1 — cos(/iAj) — cos(/cAj) + cos((/i — k)Xj] 



hk 

Now we assume (j4.2p holds. Then, as in the proof of Proposition 12. H we have 

A2 d\ A2 dXj 

^ <c . ^„ , . < 00. 



j 

This imphes J^^ X'jf{X)dX < 00. By I^LT^ . (17:22]) and the dominated convergence theorem, 
we obtain 



hm Dh,k = / A|/(A)dA. 

To prove the necessity of ()4.2p . we assume /3j(7 — X^^^i ^) < 2. Then, as in the proof of 
Proposition 12.11 we have 

A? dX 

00. (7.23) 



J 



:|A.I>i(E.=ilA.|^0' 

We let /i = A; J, and use Fatou's lemma to (j7.2ip [note the integrand is non-negative] to derive 

liminfL)/ifc> / A^/(A)dA = 00, 

where the last equality follows from (j7.23p . So lim Dfi^k does not exist and the proof is 
finished. □ 
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Proof of Theorem 14.61 

If v{t) has continuous second-order partial derivatives, then Theorem 14.21 implies that X 
has mean square partial derivatives in all directions. Let Vx{t) = {X'i{t), ■ ■ ■ and 
we show that it satisfies (|4.7p . 

For any unit vector u in M^, we can write it as ti = "^jLiUjej and "^f^iuj = 1. So 

u^^xit) = Ef=iUjX'^it)- Hence 

J^it + hu)-Xit)_^,V 



1 / ^ \^ 2 ^ / 

-^v{hu)+El "^UjXjit)] --'^Uj[EX{t + hu)Xj{t)-EX{t)Xj{t) 
^ i=i ^ i=i 

1 / ^ \2 1 ^ 

1=1 ^ 1=1 



The last equality in (j7.24p follows from (j4.5p . 

Since f (t) is an even function with v{0) = and has continuous second-order partial and 
mixed partial derivatives, then Taylor's theorem implies 

hm ^.(M = hm = -u n{0)u, (7.25) 

where fi(0) is an iV x iV matrix, with (^^(0)) = v"j{0) for i / j, and (^^(0)) .. = v-'(O). 

For the second term in the last line of (j7.24p . note that for any i, j = 1, • • • , and i ^ j, 
and any / > 0, m > 0, 



El 



X{t + lei) - X{t) X{t + mej) - X{t) 
I m 



=-^E,(^X{t + lei)X{t + mej) - X{t)X{t + mej) - X{t + lei)X{t) + X'^{t)^ ^^''^^^ 
= (^t'(^ei) + v{-mej) - v{lei - mej) 

Let / — > 0, m ^ 0, then the last term in (j7.26p goes to ^v'-j{0), where ^^^(0) is the second-order 
mixed partial derivative of u at in the i-th and j-th directions. By Theorem 14.41 we have 
E(Xj(t))^ = ^v'-{0), for j = 1, • • • , N. Hence 



. N \2 1 

E(^5^n,X^(t)j =-u^niO)u. 
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Finally for the last term in (j7.24p . we use Taylor's theorem again to derive 



1 ^ 

lim — Uiv'Ahu) = u^Vl{{))u. 



Combining this with (j7.25p and (j7.26p shows that (j7.24p goes to 0, as h ^ 0. This finishes the 
proof. □ 



Proof of Theorem 14.81 

Under (|4.2p . Corollary 14.31 ensures that the mean square partial derivative Xj{t) exists. In 
order to show that ^jit) has a continuous version, by Kolmorogov's continuity theorem or 
general Gaussian theory [cf. Adler (1981), Adler and Taylor (2007)], it is enough to show there 
exist constants C20 > and > such that 



Recall that 



Cis,t) 



Thanks to (j4.2p . we derive 

dC{s,t) 



dsj 



and 



So 



E[X'j{s)-X'j{t)]^ <C2o\s-t\^. 

(e'(^.A> _ l)(e-*<*'^) - l)f{X)dX 
cos (s — t, X) — cos (t. A) — cos (s, A) + 1 f{X)dX. 



Xj sin {s — t, X) + Xj sin (s, A) /(A) dX 



A^ cos (s -t,A)/(A) dX. 



dC{s,t) 
dsjdtj 



/ A2(l-cos(sfc-tfc)Afc)/(A)o(A 

A^(l - cos(sfc - tk)Xk)f{X)dX 
[l - cos{sk - tk)Xk)X^, 



|A|>1 



< C21 \s - + C22 ^ 
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Q+2 



dX. 



(7.27) 



E(Xj(.) - X'^it))' = E{X'^is)f + E{X'^it)f - 2E(Xj(.)Xj(t)) 

A^(l-cos {s -t,X))f{X)dX. 

The rest of the proof is similar to that of Lemma 13.21 Denote sq = i, si = (si, ^2, ■ ' ' ^^n), 
h = {si,S2,t3, ■ ■ ■ ,tM), ■■■ , SN-1 = (si, • • • ,SN-i,tN) and SN = s. Then 

N 

E(Xj(s) - X'jit))^ < iV^E(Xj(4) - X'^ih-i))^ 

k=i 



(7.28) 



Now we estimate the last integrals in (j7.28p . For simplicity of notation, we only consider 
the case when k = j [the cases oik ^ j are similar]. Denote hk = Sk — ifc and, similar to (j7.4p . 
(j7.5|) and (j7.6p . we derive 

(1 - cos(sfc - tk)Xk) Xl 
^ n f°° tu \ ^^\2J^ I dXi ■ ■ ■ dXk-idXk+i ■ ■ ■ dXN 



+ 2 1 - cos 



{hkXk))Xf^ dXk j ^ dXjg 



/■°° (l - cos(/ifcAfc))A? ^, , 

<c 2im dXk + c / Xk{l - cos{hkXk))dXk 

<C23f|/ifcp(^^-'nog-L + |/j, 
\ \tlh\ 



thanks io Hk> \. Combining this with (j7.28p proves (|7.27p . 

It follows from (|7.27p ) that the Gaussian field X'- = {X'-{t),t G R^} has a continuous 
version [which will still be denoted by X'-]. Now we define a new Gaussian random field 

X = {X{t),t£ R^} by 



X{t) = X{ti,--- ,tj.i,0,tj+i,--- ,tN)+ / XUh,--- ,tj-i,Sj,tj+i,--- ,tN)dsj. (7.29) 

Jo 

Then we can verify that X is a continuous version of X and, for every t G M^, Xj{t) = Xj{t) 
almost surely. This amounts to verify that for every t G M^, 

E{X{tf)=v{t) and E\{X{t) - X{t))' 

which can be proved by using (j7.29p . Theorem 14.41 and (j4.5p . Since the verification is elemen- 
tary, we omit the details. This proves Part (i) of Theorem 14.81 

It remains to prove Part (ii) of Theorem 14. 8[ By applying Part (i) to j = 1, we derive that 

there is a continuous version X(^) of X such that ^^^{t) is continuous. Then we apply Part 

(i) to XW with j = 2 and obtain a version X^^) of X(^) defined by 

X(^){t) = X(^){ti,0,t3,--- ,tN)+ x('\{ti,S2Mr-- ,tN)ds2. (7.30) 



Then ^^^^ (t) and ^^^^ (t) are almost surely continuous. Repeat this "updating" procedure 

for J = 3, • • • , A^, we obtain a continuous version X(^) of X such that all first-order partial 

derivatives of X(^) are continuous almost surely. Hence the sample function of X(^) is almost 
surely differentiable in the sense of (I4.10p . The proof of Theorem 14.81 is finished. 
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Proof of Proposition 16.11 

By stationarity, we'll have 

E{X{x, t) - X{y, s)f = ^X{x, t)f + ^X{y, s))^ - 2E(X(x, t)X{y, s)) 
= 2C(0,0) - 2C{x -y,t- s). 

And 

2^2 / c|x|2t 



2C(0, 0) - 2C(x, t) = 2a' - ^ ^^^^,^yr./, exp ( 

(1 + a|t|2")/5^/2 - exp 



(1 + a|t|2")/5^/2 (l + a|t|2")/37. 

cjxpT ^ (7.31) 

2 X (l+a[t|2°)/37. 



(1 + a|t|2")W2 

By using Taylor expansion, we can write (j7.31|) as 

1 + Maltpa I of |t|2"l _ 1 + ^1^1^^ o( ^l^l^"^ 

(l + alt|2'^)W2 

(l + a|t|2«)W2 

Hence we can find positive constants C24 < C25 such that 

C24(|xp^ + |tp") <2C(0,0)-2C(x,t) <C25(|x|2^ + |t|2") (7.32) 

for all (x,t) G R'^+^ with |x| and |t| small. Replace x and t in (j7.32p with x — y and t — s 
respectively, then (j6.2p follows. 

To prove (j6.3p . we make use of the fact that for any Gaussian random vector {U,V) with 
mean 0, 

(Puy - (^1/ - ^v)^) ((^f/ + ^v)^ - Pu,v) 



Var(C/|y) 



4^2, 



where p^y = E[{U - Vf], = E([/2) and ct^ = E(y2). Let U = X{x,t) and V = X{y,s), 
we derive 

,^ [C(0, 0) -C{x-y,t- s)] [CjO, 0) + Cix-y,t- s)] 

\aT{X{x,t)\X{y,s)) = ^ 

>C24\x-y\'^ + \t-s\'"). 
This proves (j6.3p . □ 

Proof of Proposition 16.21 

Eq. (|6.4p follows from Proposition 16.11 and Theorem 15. 1[ Then let's prove ()6.5p . where 
< a < 7 < 1. By Proposition 16.11 and Theorem 15. H we get 

dimGrX([0, Ij'^+i) = ^^min^^ |^|^ + d + l- A: + (l- H,)p; E 

.7 = 1 .7 = 1 
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where Hi = a, H2 = ■ ■ ■ = Hd+i = 7- Denote 



S{k) = Y,^ + d + l-k+{l-Hk)v- 
i=i 

We have 5(1) = d + I + {I - a)p, S{k) = d + ^ + {I - -f)p = S , ior 2 < k < d + 1. Also 
S?;} ^ = ^ + f Since 

p<- ^ {l-a)p< ^ 5(1) < 5 
a a 

-, d+i -, (7.33) 

^ d+^ + (l-7)p<d+- ^ ^<E^' 



and 



a a , Hi 

.7=1 J 



07 07 , Hi 

one can see that (|6.5p follows from (j7.33p and (|7.34p . 

If < 7 < a < 1, then = • • • = Iff ^ = 7 and Hd+i = a. So 5((i+ 1) = ^ + 1 + (1 - a)p, 

S{k) = d + 1 + (1 - 7)p = 5, for 1 < A: < d and Yfit\ ^ = l + i- Similarly, by comparing 

these three terms, we will get (|6.6p . □ 

Proof of Proposition 16.31 

By Proposition 16.11 and Theorem 15.21 we get that when ^ + ^ < for every x ^ W, 
y^r^{x) = a.s. And also, when ^ + ^ > Pi then for any x G M^, with positive probability 

dim(X-^(x)) = min + d + I - k -Hkp\ ■ 

^ " i<k<d+i\j^^Hj ""^j 

If < a < 7 < 1, we have Hi = a, H2 = ■ ■ ■ = H^+i = 7- Denote 



T{k)=Y,^ + d+l-k-HkP, 

then r(l) = d+l -ap, T{k) = d + ^ - 7p = T, for 2 < fc < d + 1. Since r(l) < T p < 
(16771) follows. 

IfO<7<a<l, then Hi = ■ ■ ■ = Hd = 7 and Hd+i = a. It follows that T[d + 1) = 
^ + 1 - ap and T{k) = d + I - jp = f, iov I < k < d. Since 

f<T{d+l) ^ p<-, 

7 

we derive (|6.8p . The proof is finished. □ 
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